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A planar phase space having both position and momentum non- 
commutativity is defined in a more inclusive setting than that 
considered elsewhere. The dynamics of a particle in a gravita- 
tional quantum well in this space is studied. The use of the 
WKB approximation and the virial theorem enable analytic dis- 
cussions on the effect of noncommutativity. Consistent results 
are obtained following either commutative space or noncommu- 
tative space descriptions. Comparison with recent experimental 
data with cold neutrons at Grenoble imposes an upper bound on 
the noncommutative parameter. Also, our results are compared 
with a recent numerical analysis of a similar problem. Finally, we 
provide a noncommutative version of the virial theorem for the 
case at hand. 



1 Introduction 



The subject of noncommutativity has a long history culminating in the pa- 
per by Snyder y^. The hope was that noncommutative spaces would pro- 
vide a natural background for a possible regularisation of quantum field 
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theories^ 12]. Since then different structures of noncommutative space have 
been considered [21 IH IS]- Perhaps the most widely studied case is the canon- 
ical description; i. e. when the noncommutative parameter is a constant. 
Theories defined on such a noncommutative space are based on the Weyl- 
Wigner correspondence, in which all products are replaced by the star prod- 
uct. 

A particularly interesting physical manifestation of canonical noncom- 
mutativity occurs in the context of the Landau problem-the planar motion of 
a charged particle subjected to a constant magnetic fieldjnilZ]. In refs. [HIHl 
it was shown how noncommutativity could be shifted from the coordinates 
to the momenta and vice-versa. Also, the implications of noncommutativity 
in both phase space and configuration space variables were discussed. 

In this paper we study the phenomenology of a quantum mechanical 
model with constant noncommutativity in both coordinates and momenta. 
This model is applied to the problem of a particle in the quantum well of the 
Earth's gravitational field to determine the effects of noncommutativity on 
the energy spectrum. Our analysis is totally analytical and compared with 
the numerical approach to a similar problem which has been done recently [S]. 
Finally we use the experimental results of jTUl E] to put an upper bound 
on the noncommutative parameters. We also observe the asymmetric be- 
haviour of the noncommutative parameters appearing among coordinates or 
momenta. 

This paper is divided in five sections. In the next section we define 
the noncommutative space. A general phase space transformation is given 
that maps the noncommutative space variables with the commutative coun- 
terparts. Different parametrisations are discussed which reproduce the non- 
commutative algebra quoted in the literature[71 |Hl E]. We also show that, 
contrary to recent claims [0], Planck's constant h need not be modified in 
order to simultaneously obtain noncommutativity in both coordinates and 
momenta. In section 3, the quantum gravitational well is introduced. After 
summarising the theoretical and experimental results of the energy spec- 
trum in usual commutative space, we define the corresponding Hamiltonian 
in noncommutative space. The structure of this Hamiltonian is explicitly 
obtained in both noncommutative and commutative descriptions, leading to 
completely equivalent results. In section 4 the energy spectrum is computed 
by using the WKB approximation. Compatibility with the results obtained 
from the virial theorem is shown. An upper bound on the noncommutative 
parameter is derived by comparing with the recent experimental results given 
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in pi]. Concluding remarks and discussions are given in section 5. There 
are two appendices; in the first we discuss rotational symmetry on the non- 
commutative plane while in the second we derive a noncommutative version 
of the virial theorem for our problem. 



2 Noncommutative Phase Space 

In this section we describe the noncommutative phase space on which the 
model would be defined. Both configuration and momentum space noncom- 
mutativity are considered. It is generally believed that this type of non- 
commutativity might lead to a redefinition of the Planck constant. However 
this is not mandatory, as our construction shows. 

Consider a two dimensional space, where the position and momentum 
operators satisfy the standard Heisenberg algebra 

[xj , Xj] 

h,Pj] = (1) 
[xi,pj] = ihSij. 

Note that this algebra is invariant under the following symmetry transfor- 
mation 

■^i ^ Pi 

Pi Xi (2) 
i —>■ —i. 



It is then possible to develop a coordinate space representation where the 

dxi 



coordinates Xi are diagonal and pi = —ih^^- or, alternatively, a momentum 



space description where the momenta pi are diagonal while Xi = ih^r-. 
Now we take the following general phase space transformation 

yi = Xi + aieijPj + a2eijXj (3) 
qi=Pi + PieijXj + (52^ijPj (4) 

where a, /3 are arbitrary constants. Here we enforce a symmetry leading to 
Ui Qi and Qi — > Hi under the transformation Clearly this is possible if 
we introduce the following transformation 

ai f3i (5) 
A ^ «i- (6) 
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Thus a symmetry transformation, analogous to (j21), in the modified y — q 
plane is given by 



qi Vi 

ai Pi (7) 

A 

i — i> —i. 



Using we can show that the new coordinates y and momenta q satisfy 
the algebra 

[yi^yj] = -2ihaieij (8) 

[qi,qj] = 2ih(3ieij (9) 

[yi, qj] = ^^(1 + "2/^2 - ttiA)5ii + «^(«2 - P2)(^ij- (10) 

Under the symmetry transformation ((7j) the above algebra is invariant. 

So far we did not associate any specific values to the coefficients a and 
/3's. Now if we set 

9 



2h 
JL 

a2 = (32 = 



A 



we obtain 



[yi,yj] = i^f^ij 

[qi, qj] = irjeij (11) 

[yi, qj] = ^(1 + = i^effSij 

which reproduces the noncommutative structure given in 0. The term ^ 
is interpreted [H] as the correction to the Planck constant. This correction is 
claimed as a consequence of noncommuting phase space algebra. As men- 
tioned in [121 it is possible to scale 9 and rj in a proper way so that heff = h 
but in that case the first two relations of eq. (fTTj) take rather complicated 
forms and the noncommutative parameters cannot be expressed simply as 9 
or T]. However it has been shown by us [HIH], using a different approach that 
both 9 and r] noncommutativity can be retained without altering the Planck 
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constant. This is also achieved in the present context by taking the following 
values of a and j3 

"1 = —k 

= (12) 
which yields the noncommutative algebra 

[yi,yj] = 

[qi,qj] = irjeij (13) 
[yi,qj] = ihSij 

so that the Planck constant is not modified. Physical applications of this 
type of noncommutative algebra may be found in [7| |H] . 

It should however be mentioned that there exist studies ^2] which show 
that the effect of a modified Planck constant was negligible with respect to 
the magnitude of the effect of interest. In this sense the modification (or 
otherwise) of the Planck constant is more of academic interest. 

The inverse phase space transformation is given by 



where 



Xi = Ayi + Beijyj + Cqi + Deijqj 
Pi = Eyi + Feijy-j + Aqi + Btijqj 



(14) 



2{h? - Or]) ' 2(^2 - 9rj) ^ ' 



E 



Observe that 9 and i] must have different signs so that the various coeffi- 
cients are real and well defined which guarantees the hermeticity of physical 
operators x,p and y, q. 

3 Gravitational Well 

The problem is first discussed in usual commutative space. We consider a 
two dimensional plane where a particle of mass m is subjected to the Earth's 
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gravitational field in one direction; the vertical taken to be described by the 
coordinate Xi. We assume that the gravitational acceleration g is constant 
near the surface of the earth. The commutative Hamiltonian is given by 

H = ^{pI + pI) + mgxi. (16) 

Since the particle is free in the X2 direction, its energy spectrum is continuous 
in that direction and the wave function can be written as 

^(X2) = j g{k)e'^''^dk. (17) 

In the other direction the wave function is the well known Airy function 0((^) 
with appropriate normalization|14j, 

Mx{) = A^m ; e = ' (^1 - — ). (18) 

\ h'^ J mg 
The zeroes of the Airy function, jSn give the energy eigenvalues 

K = -(^^y/5n; rz = l, 2, 3... (19) 

Below the classical turning point a;„ = ^ the wave function oscillates and 
above Xn it decays exponentially. This was observed experimentally by 
Nesvizhevsky et a/.[10j. They form the "well" by placing a horizontal re- 
flecting mirror in the Earth's gravitational field. The neutron was used 
as the quantum particle since it is chargeless and has a longer life time 
(r ~ 885.7s) |15J. By placing an absorber above the mirror they allow a cold 
neutron beam to fiow with a horizontal velocity V2 = 6.5ms~^ between the 
mirror and the absorber. Then they measure the number of transmitted 
neutrons as a function of absorber height: this was shown to be a step like 
function which manifests the quantum nature of the problem. 

For theoretical computations, a more transparent and simple solution is 
obtained in the WKB approximation. The energy eigenvalues for the linear 
gravitational potential are given by. 

En = (^— [7r%(n--)]2j (20) 
= ang-^ ; n = 1, 2, 3... (21) 
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where 

'9m, ^, 1 



«n = \^VKr^'f\) ■ (22) 

A summary of both theoretical and experimental results is given. Taking the 
values of constants as 

h = ^(Planck constant) = 10.59 x 10^^^ Js (23) 

g = gravitational acceleration = 9.81 ms~^ (24) 

m = mass of neutron = 167.32 x 10^^^ Kg (25) 

the first two energy levels found from (|2Up are^, 

El = 1.392 peV = 2.23 x lO^^^J (26) 

E2 = 2.447 peV = 3.92 x lO^^^J. (27) 

From El and E2 the classical turning points are calculated to be 

El , , 

xi = — = 13.59/im 28) 

mg 

xx = — = 23.88/im. (29) 
rag 

These are in reasonable agreement with the experimental results [TT] 

P = 12.2 ± 1.8(syst) ± 0.7(stat) (^um) (30) 

x^^ = 21.6 ± 2.2(syst) ± 0.7(stat) (/im). (31) 

Error bars for the above mentioned energy levels are 

AEf P = 6.55 X 10-^2 J = 0.41 peV, (32) 

A£;f P = 8.68 X 10"^^ J = 0.54 peV. (33) 



^The error is ^ 1% compared to the results derived from H19|) 
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3.1 Noncommutative Space Description 

The problem is next formulated on a noncommutative space. The analogue 
of the Hamiltonian ()16|) in noncommutative space is defined as 

H = ^{ql + ql) + mgyi (34) 

where the variables satisfy the algebra pHj) . To find the spectrum, two 
approaches are possible. One can directly work in the noncommutative space 
variables or use the phase space transformations to reduce the problem on 
the usual commutative space. We first discuss the second approach. Using 
the maps ()3|4|) together with the parametrisation p2j). we find, 



ff , 2 2\ VV~^V / 2 2\ 

(Xi + X2) + — r:r{^iPi + Pi^i) - o ^iPl+P2)- (35) 



8mh? 8mh? 8mh? 

Defining a new constant 

and a new variable 

P2=p2- m^g-f (37) 
we can write the above Hamiltonian in the form 

+ ~5 ^(^l^"! + ^2P2 + PlXl + P2a;2) 

8mh'^ 



+mg{{l - -)x, + -^{1 + -).,} - ^p^^. (38) 

Since the difference between p2 and p2 is just a constant, they satisfy the 
same commutation relations. The eigenvalues of p2 are translated by an 
equal amount vis a vis those for p2 and hence these are not distinguished. 
Also neglecting the additive constant in the Hamiltonian (jHHj) we get 



This Hamiltonian is put in a more transparent form by first defining a mod- 
ified gravitational acceleration g in the following way 

mg cose = mg{l — — ) (40) 

mg sme = mg (1 + — ). (41) 

The tilting angle with the xi axis is given by 



2h \2h-r]-f 

while, 



1 



Since the product Orj is negative, g is always positive definite. Now we rotate 
in the Xi — Xi plane by an angle e, so that the coordinate of a point in the 
rotated frame is given by 



x\ = cose xi + sine X2 
x'2 = cose X2 — sine Xi. 

Correspondingly, the momenta are transformed : 

p[ = cose pi + sine p2 
p'2 = cose p2 — sine pi. 



(44) 



(45) 



Using and it is easy to show that 

p'i+P2=pI+pI (46) 

x[p2 - x'^p'i = X1P2 - X2P1 (47) 

xf + x'i = xl + xl (48) 

x'iP'i + p'ix'i = XiPi + PiXi. (49) 

Therefore in the rotated frame the noncommutative Hamiltonian is given by, 



x[p-+pWi)-^iPi+P2)- (50) 
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The primed and unprimed variables satisfy the same algebra; henceforth 
the primes are all dropped. Then we can identify the first three terms of the 
Hamiltonian (jSUj) exactly as the commutative Hamiltonian given in p6|) . This 
should be considered as the unperturbed Hamiltonian. The term ^^^ijPiXj 
is effectively a Landau problem like term, where a magnetic field is present 
perpendicular to the Xi — X2 plane. The term g^^(a;f + x^) is practically an 
oscillating potential. 

Since the noncommutative effects are rather small we first confine to 
the leading order approximation in 9 and r]. Moreover ()43|) shows that in the 
leading order, 

~g = g[l + 0{er])]. (51) 
Hence the Hamiltonian (j5Uj) in the first order approximation is given by, 

H = ^(Pi + P2) + 1^9x1 - ^^(^^P2 - (52) 
= Ho - ^^{xiP2 - X2P1) (53) 

where Hq is nothing but the commutative Hamiltonian already given in eq. 
fll6|) . The energy spectrum pertaining to this Hamiltonian will be computed 
in section 4. 



3.2 Alternative Formulation 

Here we analyse the structure of the Hamiltonian directly in terms of non- 
commuting space variables. In the leading order approximation, shows 
that only the noncommutativity in the momentum variables is relevant. Thus 
to simplify matters, consider 9 = from the very beginning. It is then possi- 
ble to provide a representation of the operators in the coordinate basis, such 
that the coordinates yi are diagonal while the momenta Qi are given by the 
usual differential operator plus an extra piece, 

Vi Vi (54) 

-'^w^ + i'^^'^- ^^^^ 

It is simple to verify that this is a valid representation of the algebra (fT^ 
with ^ = 0. In this basis the Hamiltonian ()34|1 takes the differential form 
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It is now possible to compare this with the expression given in (j5(J|) . Since 
(|5(J|) involves commutative space variables, either the coordinate or momen- 
tum representation may be used. Specifically, in the standard coordinate 
representation one may replace x[ yi and p'^ — >• —ih-^. Then for ^ = 
the result is reproduced. Obviously in the leading order approximation 
this leads to 

It is straightforward to repeat the analysis for the case = 0. Ex- 
pression ()50|) reveals that there are no corrections and only the unperturbed 
Hamiltonian remains. Coming back to the noncommutative variable ap- 
proach, the algebra ()13|) implies that, for = 0, a momentum space repre- 
sentation can be constructed. This is given by, 

qi Qi (57) 
d 9 



The Hamiltonian (1341) is then written as 



= ^(«r + «l) + -|:*). (59) 

This is next compared with ()50|) for 77 = 0. Apparently it seems that there 
is a mismatch. Using the standard momentum representation in ()50p so that 
p[ Qi and x[ results do not agree. However, as has been 

stressed even at the basic level by Dirac^H], both coordinate and momen- 
tum representations admit a generalisation. This generalisation is, however, 
trivial since it can be absorbed in the phase of the wave function. In this 
problem we can set p'^ Qi and x[ ^^^~^Q'2 because the algebra remains 
unaffected. With this choice (jMIj) agrees with (j^ . Thus the commutative 
and noncommutative space descriptions are equivalent up to an irrelevant 
over all phase. 

We conclude this section by observing that a direct comparison between 
commutative and noncommutative space descriptions is possible only if either 
^ = or ?7 = 0. Then as can be seen from ()13p . it is feasible to formulate the 
coordinate representation or the momentum representation, respectively. If 
both 6 and 77 are nonvanishing then it is somewhat problematic to construct 
an explicit representation directly in the noncommutative space. In that case 
the variable change approach leading to is necessary. 
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4 Bounds on Noncommutative Parameters 



Here the energy spectrum is computed and therefrom bounds on the non- 
commutative parameters are determined. Consider the Hamiltonian ()53p in 
the first order approximation. Now the term proportional to r] in the above 
Hamiltonian can be treated perturbatively. The unperturbed part Hq is 
known to be exactly solvable in terms of Airy functions jl4j. Furthermore, 
using the property that Airy function (or any bound state wavefunction vis 
a vis motion in the direction xi) is real, it is easily seen that 

r+oo Q 
<Pl>n= dXiij*^{-th—iJn) = 0. (60) 

Jo 

This can also be understood physically from the fact that, for a bound state 
system, the average current flow in a particular direction is zero. So effec- 
tively the Hamiltonian turns out to be 

In this way we see that, in the leading order, the noncommutative corrections 
are entirely encoded in the term 

H, = -^..P. (62) 

Such a correction term was also found in the approach of |9j. The energy 
corrections due to this term were calculated in ^ from the standard pertur- 
bation formula 

r+co 

AEn = / dxi'ilj!^{xi)Hiiljn{xi). (63) 
Jo 

where i^nixi) is given in (fTHj) . This integral was evaluated by using numerical 
methods to obtain the following corrections for the lowest two energy levels 

|A^i| = 2.83 X lO^^r^ (J) (64) 
|A^2| = 4.94 X lO^^T] (J). (65) 

Comparing the above result with (fH^ and (j^Sl) the upper bound obtained 
on the noncommutative parameter i] was found to be [HI, 

\r]\ < 2.32 X lO"*^^ kgV^s'^ (n = 1) (66) 
\r]\ < 1.76 X 10"^^ kg^m^s"^ (n = 2). (67) 
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In our approach we avoid the numerical analysis totally and follow the 
semi classical WKB method which works extremely well for a linear potential. 
We write the complete Hamiltonian in the form 

X 77 

^ " 2^^^'+^2) +"^(^^- ^^^^^'2)3:1 (68) 
= T^bi + pI) + ^a'xi (69) 

where g' = fl'~2r^^'2- Since in the X2 direction the particle is free, p2 is a con- 
stant of motion. In the experiment painstakingly performed by Nesvizhevsky 
et al. |1Q^ the expectation value of p2 was 

< >= 10.91 X 10"^^ Kg m s"\ (70) 

Now we can use (12111 to write the corrected energy values of the Hamiltonian 

as, 

where En corresponds to the unperturbed energy and AEn is the correc- 
tion. It is possible to find an analytic expression for Ai?„ from ()71|) by an 
expansion, 

K + AK = a„<75(l--^ <P2 >)i (72) 

Retaining the leading r^-order term we find, 

= -TT^ < P2 > En. (73) 

As expected the same functional form of the result follows through the use 
of the virial theorem [T7]. 

For a linear potential the virial theorem implies < T >= ^ < y > 
where T and V are kinetic and potential energies, respectively. Then the 
total energy E is given by E =< T > + < V >= ^ < V > . In our 
problem V = mgxi, so that < xi >= Now the perturbation term is 

AE = < P2 >< Xi >. Hence we find 

AE = -TT^ <P2> (74) 
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which reproduces the structure (j73|) . Taking the values of Ei and E2 from 
EZD and <P2> from ((701) we get on using ((73), 

|A^i| = 2.79 X lO^^r; (J) (75) 
|A£^2| = 4.90 X 10^^77 (J). (76) 

Finally, using the experimental input from (jH^ IHHj) leads to the following 
upper bounds on rj; 

\r]\ < 2.35 X 10"^^ kg^m^s-^ (n = 1) (77) 
\ri\ < 1.77 X 10"^^ kgV^s~2 (n = 2) (78) 

The energy corrections ITBj) are in excellent agreement with the numer- 
ical results obtained by perturbing about the exact Airy function 
solutions. The same naturally holds true for the upper bound on rj. 

5 Conclusions 

We have discussed a model of a particle in the quantum well of the Earth's 
gravitational field and a perfectly reflecting horizontal plane beneath, de- 
fined in a space with noncommuting coordinates and momenta. The energy 
spectrum in this model was computed analytically by exploiting the WKB 
approximation. Comparison with the experimental findings of |im ITT] placed 
an upper bound on the rj- noncommutativity parameter appearing in the al- 
gebra of momenta. 

Our results were also in excellent agreement with the numerical calcu- 
lations of in] where a similar model was considered. However to put things 
in a proper perspective, we emphasise that there is a conceptual difference 
between our model (j34p and that of |9., . In our treatment the model is defined 
directly in the noncommutative space. It is then possible to analyse it either 
in terms of noncommuting variables (section 3.2) or, using the phase space 
transformations, in terms of commuting variables (section 3.1). The results 
agree in either formulation. In ref. P], on the contrary, the model was de- 
fined by taking the usual Hamiltonian p6|) and exploiting the inverse phase 
space transformation (jl4|) to express it in terms of noncommuting variables. 
A perturbative expansion of the Hamiltonian in i] and 9 is next carried out 
whereas we get a closed form expansion (I^Uj) which truncates at the second 
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order in the noncommutative parameters. For the leading order term, how- 
ever, the two expressions match, upto a sign of r]. This permitted us to make 
the necessary comparison. 

An essential ingredient of this analysis was to exploit rotational symme- 
try to simplify the Hamiltonian. Normally this symmetry would be violated 
for constant noncommutativity. However, as shown in appendix 1 for the 
particular case of two dimensions, it holds. 

Lastly, in appendix 2, we derive a noncommutative version of the virial 
theorem for the linear gravitational potential considered here. 

Appendix 1 

Here we discuss the validity of rotational symmetry in two dimensions. To 
see this take the following general structure for arbitrary n-dimensions, 

[yi,yj] = idij, ij = 1,2, ...n (79) 

Under rotations, 

% = i^ijUj ; = -^ji- (80) 

Then taking the variation on both sides of ()79|) . 

[Syuyj] + [yuSy,]=0 (81) 

which implies, 

^ikOkj + ^jk^ik = 0. (82) 

This is not true in general. For d = 2, we may write ooij — udij, Oij — ^^ij so 
that the above condition simplifies to, 

^^{^ik^kj + ^jk^ik) = (83) 

which holds. This shows that rotational symmetry is valid in d = 2 noncom- 
mutative space. For higher dimensions deformed transformations have to be 
appropriately defined which restore this svmmetrv |18t IT^ I2(jj. 
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Appendix 2 



A noncommutative version of the virial theorem in the present context is 
presented here. For a hnear potential, the usual virial theorem states that 
the average kinetic energy of a particle is half of the average potential energy. 
In a noncommutative space this theorem gets modified. Using the Heisenberg 
equation of motion, for any arbitrary state we can write 

^ < Viqi >= ^ < [ViQi^ H] > . (84) 

where i/i, Qi are the noncommutative space variables. For our case we take 
the Hamiltonian H as given in ()34|) . For stationary states the left hand side 
of the above equation is zero. Using the explicit commutation relations given 
in (fT!^ . the eq. turns out to be 

77 6 

2 < — > - < mgyi >= < eijyiQj > +-mg < q2 > ■ (85) 

2m mh h 

This is the modified virial theorem where the corrections are found to be first 
order in 9 and rj. In the limit 6,1] —>■ the usual virial theorem is recovered. 
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